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Abstract 

Fill an n X n matrix with independent complex Gaussians of variance 1/n. As 
n — > oo, the eigenvalues {zj.} converge to a sum of an //^ -noise on the unit disk 
and an independent H^^^-noise on the unit circle. More precisely, for functions 
of suitable growth, the distribution of Ylk=i(f i^k) — ^/(zfc)) converges to that of a 
mean-zero Gaussian with variance given by the sum of the squares of the disk and 
the circle H^^^ norms of /. As a consequence, with p„ the characteristic polynomial, it 
is found that log \pn\ — Elog \pn\ tends to the planar Gaussian free field conditioned to 
be harmonic outside the unit disk. Further, for polynomial test functions /, we prove 
that the limiting covariance structure is universal for a class of models including Haar 
distributed unitary matrices. 

1 Introduction 

Consider the Ginibre ensemble, that is the n x n random matrix in which all entries are 
independent complex Gaussians of mean zero and variance 1/n. The eigenvalues zi, Z2, . . . , Zn 
form a point process, a realization of which is depicted on the left of Figure ^ A first glance 



Figure 1: Ginibre eigenvalues and uniform points in U, n = 150. 
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Figure 2: The running eigenvalue sums z^, m = 1, . . . , 150, in which the Zk are ordered 
by their arguments, and the same for the independent points. 

at this picture suggests that the eigenvalues are uniformly distributed in the unit disk U. 

Indeed, Bai has shown that, with probability one, the empirical measure converges to the 

uniform distribution on U. This is the circular law of the title. 

Compared to n points dropped independently and uniformly on U (see the right of Figure 
the Ginibre points are clearly more regular. For example, their sum, the matrix trace, 

is complex Gaussian with variance 1, while for the independent points the variance is |n 

(see Figure 121 for comparison). Our main theorem addresses this phenomenon in greater 

generality, rigorously verifying a prediction of Forrester [7j. 

Theorem 1 (Noise limit for eigenvalues). Let / : C ^ M possess continuous partial 
derivatives in a neighborhood of U, and grow at most exponentially at infinity. Then, as 
n oo, the distribution of the random variable '^^=1 ifi^k) — Ef{zkj) converges to a normal 
with variance crj + aj, where 

is the squared Dirichlet (H^) norm of f on the unit disk, and 

~^] = \\\f\\mn^3V) = \Y.\^\\m' 

is the H^/'^-norm on the unit circle dV, where f{k) = ^ J^^ f(6)e~^^^d6 is the k-th Fourier 
coefficient of f restricted to\z\ = l. 

We mention that the centralizer Ylik=i ^fi^k) is easily replaced by n x ^ J^^ f{z)d'^z with 
no change to the outcome. More interesting is the covariance structure that appears above, 
being composed of two independent terms: one for the boundary and one for the bulk. An 



interpretation of this limiting noise is provided upon considering the characteristic polyno- 
mial Pn{z) = 11^=1 (-^ ~ ^k) which has fluctuations described in terms of the Gaussian free 



The planar Gaussian free field (GFF) is a model which has received considerable recent 
attention as the scaling limit of uniformly random (discrete) ^ M surfaces, though it 
apparently has not previously been connected with any matrix model. While the realizations 
of the GFF fluctuate too wildly to allow it to be defined pointwise, it may be defined as 
a random distribution h as follows: for all suitable test functions /, the random variable 
(/, is Gaussian with variance = /c I V/P '^^-2^- more classical language 

this is the Gaussian Hilbert space for H^{C), identifying the GFF as the most natural 2- 
dimensional analogue of Brownian motion. 

As the Ginibre eigenvalues accumulate in the unit disk, log|p„(z)| tends to be harmonic 
away from the unit disk as n — > oo. Thus the following version of the GFF appears naturally 
in the limit. 



Corollary 2 (The Gaussian free field limit). For the Ginibre characteristic polynomial 



Then hn converges weakly without normalization to h* , the planar Gaussian free field condi- 
tioned to he harmonic outside the disk. More precisely, for functions f as above we have 



in distribution, and the same holds for the joint distribution for the integrals against multiple 
test functions fi, . . . , f^. 

Next, by the Cramer- Wold device, a version of Theorem ^ holds for complex valued test 
functions / of like regularity and growth. In this case the noise term is replaced by the 

norm of the d = {l/2){dj: + dy) derivative of /. Then, if it is further assumed that / is 
complex analytic in a neighborhood of U, the covariance structure simplifies significantly: 



where A/c denotes complex Gaussian distribution with given mean and variance. This should 
be compared with the result of Rider and Silverstein ^H]; where the same central limit 



field. 



hn{z) = log \Pn{z) \ -Elog \Pn{z)\. 




J2ifi^k)-Efiz,)) ^ Arc(o,-||/'||i.(^)) =Arc(o,||/||^v2(af/)) 



k=l 
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theorem is proved for random matrices with more general than A/c(0, 1) entries, but where 
the test functions need to be analytic throughout \z\ < 4. 

Perhaps surprisingly, exactly the same limiting covariance structure arises from eigenval- 
ues of Haar distributed unitary matrices (Diaconis and Evans, [H]). While the difference of 
the two limiting noises seems very clear (one lives on the disk and one on the circle), analytic 
functions do not see this. We show that this is not a coincidence by proving a universality 
result. 

The important point is that both the Ginibre eigenvalues and those of Haar U (n) comprise 
symmetric polynomial projection (SPP) processes in C. An SPP process is a determinantal 
process defined by a pair {fi,n), where /i is a rotationally invariant probability measure on 
C and n is the number of points, as follows. For k = n the joint intensity p{zi, . . . , Zk) with 
respect to fi'' is given by det{K{zi, Zj)i<ij<k), for K the integral projection kernel to the 
subspace of polynomials of degree less than n in L^(C,/i). In fact, as soon as this holds for 
k = n, it holds for any k, see for example [TTj . 

Theorem 3 (Universal noise limit for analytic polynomial test functions). Consider 
a sequence of SPP processes defined by {fin,n). Assume that for all integers m we have 



in distribution. 

The proof of Theorem |21 relies on the close connection between SPP processes and Schur 
functions. For unitary matrices, this connection is a consequence of Frobenius duality in 
representation theory (see 0); in the last part of the paper we show that they are useful 
in a more general setting. For the Ginibre ensemble /i„ = A/c(0, 1/n), while for unitary 
eigenvalues /i„ is uniform measure on the unit circle. When /x„ is uniform measure in the 
unit disk we get the (truncated) Bergman ensemble, for which the n oo limit is a discrete 
point process in U that agrees with the zeros of the random power series with i.i.d. complex 
Gaussian coefficients, see jTTj. Each of these ensembles satisfies the moment condition laid 
out in the statement, and in each case the result actually holds for / analytic. 

Methods. Central limit theorems of the above type have been studied in the random 
matrix theory community for matrix ensembles Hermitian ensembles, e.g., |13| lUl 131 IT|. and 



M{n, 2n + 2m) /M{n, 2n) 1, 



where M{n, k) 




as n — > cx). Then, for f any polynomial in z, 




4 



ensembles connected with the classical compact groups, e.g., [121 123 El • In all of these 
cases, the eigenvalues form a one- dimensional point process either on the real line or on the 
unit circle. 

The situation for two-dimensional eigenvalue processes is more complicated; indeed, it 
seems that this paper is the first to prove a central limit theorem in two dimensions for 
general (not necessarily analytic) test functions. For the one-dimensional case, variants 
Wigner's method of moments argument has been refined and extended to give central limit 
theorems in several settings. However, this argument does not work for non-analytic test 
functions, as such functions are not expressible as a trace of a matrix polynomial. While 
analyticity in one dimension is essentially a smoothness assumption, in two- dimensions it is 
much more restrictive. This makes the proof of the general circular law |2] more difficult 
than that of the semicircle law, and explains why the corresponding general central limit 
theorem is still open (again, see JHj for analytic test functions). 

Our proof is based on the determinantal structure of the point process and a cumulant 
formula for determinantal processes introduced by Costin-Lebowitz and generalized by 
Soshnikov [221 12H1 121| • In the present setting, it allows us to establish a connection between 
joint cumulants of monomial test functions in z and z and a class of combinatorial objects 
which we call rotary flows. 

When the associated matrix integrals can be evaluated by separate methods, then such 
a connection can be used to count intractable combinatorial objects or evaluate difficult 
combinatorial sums arising, for example, in statistical physics. Here the cumulants become 
combinatorial sums over rotary flows. Fortunately, we are able to establish (Section ^ 
conditions under which these combinatorial sums vanish as n oo, without explicitly 
computing them for any finite n case. As it turn out, this is sufficient to conclude the 
asymptotic joint normality of monomial statistics. This result is then lifted from polynomial 
to general functionals by a concentration estimate, yielding Theorem ^ The sum of these 
steps occupies Sections jH through [71 Theorem 3 is proved in Section [HI while, after some 
preliminaries in the next section, the connection to the GFF is detailed in Section [21 

The arguments introduced here work in the general rotational invariant setting of SPP 
processes. However, beyond analytic test functions, the structure of the limit depends inti- 
mately on the process in question. In this paper we concentrate on the Ginibre ensemble. 
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2 Facts about the Ginibre ensemble 



The non-Hermitian matrix ensemble of independent complex Gaussian entries is named for 
Ginibre who discovered [S] that the joint density of eigenvalues, zi through Zn, is given by 

dHn{Zl) ■ ■■d^n{Zn) ^ JJ \Zj - ^fe^. (1) 

Here (i/i„ = ^e~'^'^'^(i^z is complex Gaussian measure of variance 1/n, and Zn < oo is a 
normalizing constant. 

A description equivalent to is to say that the eigenvalues of the Ginibre ensemble 
make up the determinantal process tied to the projection onto the subspace of polynomials 
of degree less than n in L^(C, the projection kernel with respect to fin being 

n-l 



Kn{z,w) = ^ 



That is, all finite dimensional correlation functions, or joint intensities, are expressed in 
terms of determinants of the kernel Kn'- for Bi, . . . , any mutually disjoint family of Borel 
subsets of C, 



k 

Ef 



TT#{fc:Zfce5£} = / det(Kn{zi,Zj)) d^n{zi) ■ ■ ■ d^n{zk). (2) 

One then refers to the process {Kn,Hn)- Further background information on determinantal 
processes may be found in ^1] and [22] • Immediate from either (Q) or Q is that the 
eigenvalues repel each other; Figure 1 provides a somewhat dramatic snapshot of this fact. 

As the dimension tends to infinity, the reader is invited to check by a simple calculation 
that the mean density of Ginibre eigenvalues. 



1 " 1 

\-y^.^Zkiz) = -Kn{z,z)dfin{z), 

In ^-^ J n 



E 

In 

k=l 

tends (weakly) to the uniform measure on the the unit disk U. This is a particularly simple 
instance of the circular law about which we describe fluctuations. 



3 The Gaussian free field 

For a domain D C C, let = H^B) denote the Hilbert space completion of smooth func- 
tions with compact support in D with respect to the gradient inner product {f,g)m{o) = 
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(/) ~ lo'^f ' '^g<Pz. Note that the if^-norm is invariant under conformal transfor- 
mations and complex conjugation. 

The Gaussian free field (GFF) on D is defined as the Gaussian Hilbert space connected 
with this norm. It may also be interpreted as the random distribution h such that for all 
/ G H^{p>) the random variable (/, is centered normal with 

Var ((/,/i)hi(d)) = \\f\?m{py 

For more detailed information concerning the GFF, the paper ^H] is recommended. 

In contrast, the if^-noise on D may be defined as the Gaussian random distribution h 
for which (/, h)L2 is centered normal and 

Var((/,/i)z..) = ||/||^i(„). 

The two are related by noting that 

(/,/i)Hi = (-A/,/i)i2 = (/,-A/i)^2. 

That is, formally there is the identity in law h = —Ah. 

Planar GFF harmonic outside U and Corollary 2 

The space if^(C) decomposes into three orthogonal parts (see [IH]): H^{C) = AiQ)AoQ)Ah. 
Aj is the closure (in H^) of smooth functions supported in a compact set of the open unit 
disk. Ao is the similar closure of functions supported in a compact set in the complement 
of the closed unit disk. Last, Ah is the subspace spanned by continuous functions which are 
harmonic on the complement of the unit circle. 

We can now define the GFF h on C conditioned to be harmonic exterior to U as the 
projection Vih of the standard planar GFF to Aj © Ah in the distributional sense: for 
/ G H^, set {f,ViHh) = {ViHf,h). Note that for / G and continuous, Vnf is the 
harmonic extension of the values of / on the unit circle to the whole plane. That is, the 
resulting function is harmonic at infinity, or (P/^/)(l/z) is harmonic near the origin. 

The proof of Corollary 2 requires one preliminary observation. 

Lemma 4. For Vih the projection unto the Aj © Ah it holds that 

\YPiHf\\]{i{c) = \\f\\m{v) + 2ll'^^''/II^^Mc)- 
For f continuous in the neighborhood of dU the second term satisfies 

\\\T^Hf\\]mc) = II^h/IIhhu) = ^ll/ll/fi/2(au)- (3) 
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Proof. If / G H^{C), there exists a function g G Aq so that s = f + g is symmetric with 
respect to inversion, i.e., s{z) = s{l/z). With || ■ || denoting H^{C) unless specified otherwise, 
we write 

\\sf = WVjsf + \\Vosf + WVHsf. 

By symmetry and the conformal and conjugation invariance of H^, we have that HP/sH = 
\\Vos\\, and also = INI||'i(c\u) = PIIV2- Hence, 

WPisf = ||s||^i(u) - ^WPusf, 



and we conclude that 



1„_ „o „ „„o 1 



since g = s — f ^ Aq- This proves the first claim. 

The first equality in (jS)) follows from symmetry and conformal (inversion) and conjugation 
invariance. Note that all norms in depend only on the values of / on dU, so we may 
assume that / is harmonic in U, and drop the projection. Harmonic functions in U are 
spanned by the real and imaginary parts of z'', so it suffices to check that for two elements 
gi, g2 of this spanning set 

which is a simple computation. □ 
Proof of Corollary\^ By Lemma ID the limiting variance in Theorem can be written as 



Now recall Pn{z) = Y[k=ii^ ~ ^k), the Ginibre characteristic polynomial, and set 

f^n{z) = 7^ log \p„{z) \ - T^Elog \Pn{z)\. 
zn ATT 

Since 2^Alog|z| = the delta measure at the origin in distribution, hn{z) = Ahn{z) equals 
the centered counting measure of eigenvalues in the same sense. Thus, by Theorem 1, we 
have that for any / which is once differentiable in a neighborhood of U, the random variable 
(/;^n)L2(u) is asymptotically normal with variance ^||'^/h/||//i(c) as n — oo. In this way 
we can identify the distributional limit of hn{z) as the planar Gaussian free field, conditioned 
to be harmonic outside the unit disk. □ 
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4 Cumulants, polynomial statistics and rotary flows 



For any real- valued random variable X the cumulants, Ck{X), k = 1,2,..., are defined by 
the expansion 

logE[e^*^] = 5^i|^C,(X), (4) 

k=l 

and carry information in the same manner as the moments of X. Important here is the fact 
that the variable X is Gaussian if and only if Ck{X) = for all k > 2. 

Now, as it seems to have been first observed by Costin and Lebowitz in ^ , the cumulants 
of linear statistics in (any) determinantal point process have a particularly nice form. Staying 
in C, choose a kernel K and measure fi so that K defines a self-adjoint integral operator /C 
on -L^(C, d^). Then, if K. is locally trace class with all eigenvalues in [0, 1], {K, fi) determines 
a determinantal point process (see ^T], ^H] or [221 )• That is, there is a process with all 
A;-point intensities satisfying the identity (j21)- With X{g) = J^di^k) and {zk} the points of 
{K,n), the k-th cumulant is written, 

^ — ^ m ^ — ^ k-il ■ ■ ■ km'- 

fci >l,...fcm>l 



X 



« m 

/ \T\(9{ze))''')K{zi,Z2)K{z2,Z3) ■ ■ ■ K{z^, zi)dfi{zi) ■ ■■d^{z^). 



Along with the above has been put to important use in 

For what we do here it will be convenient to cast the cumulants somewhat differently. 
First let [k] = {1, . . . ,k}. For a function a : [k] ^ [m], and / G Q'', where Q is an algebra 
real-valued functions, define af G Q"^ by 

(^A(^)= n /^w- 

i: cr{i)=j 

Fix a functional : G"^ ^ for each 1 < m < k. For f & define 
where the sum is over all surjections (denoted by ^). Further define 

^ f i\m— 1 

T,($, /) = V Tfc,„($„, /). (6) 

^-^ m 

m=X 
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Then, for the determinantal process {K,fi), we set 

*m(/l, • • • , /m) = y fl{zi) ■ ■ ■ fm{Zm)K{zi, Z2)K{z2, Z3) ■ ■ ■ K{Zm, Zl) djli^Zi) ■ ■ ■ dn^Zm), 

and replace the formula (0) for the fc-th cumulant with the equivalent, 



Ck{X{g)) = T,{^,{g,g,...,g)). 



(7) 



This new expression, as a sum over surjections rather than over ordered partitions, lends 
itself to the combinatorial approach we take. In the following, it is explained how the 
cumulant formulas, (jHI) through ((7j), simplify in the case that the determinatal process in 
question has a rotation invariant reference measure /z in C and the linear statistic considered 
is a (weighted) polynomial in z. The next section (Sectional) establishes when the limiting 
combinatorial sums vanish. Finally, these facts are put together to compute the asymptotics 
of the cumulants of polynomial statistics in the Ginibre ensemble, establishing the CLT in 
that case. 

Rotation invar iance and rotary flows 

The following identities hinge on the rotation invariance of Ginibre ensemble. To drive this 
point home, we set things up for general determinantal point processes which share this 



Begin with a radially symmetric reference measure (i in the complex plane normalized so 
that /i(C) = 1, and define the kernel. 



where it is assumed that G [0, 1], \z\'^^djj, = c^'^ < oo, while L = oo is allowed. Again, 
the results of [12] or [22] (see also JT]) imply that {K, fi) define appropriate joint intensities. 



Now return to the formula ((Tj) for the fc-th cumulant of the linear functional X{g) = 
T.9{zk): 



Since T($, (/i, . . . , fk)) is a /c-linear symmetric functional, it follows that if (7 is a polynomial, 
g = 'Y^ akfki then Ck{X{g)) is a sum of terms of the form T($, (/i, . . . , fk)) in the monomials 
fk- Our goal is to understand the conditions which will make Tfc($, (/i . . . , f^)) vanish (or 
vanish in the limit of some parameter, soon to be the dimension for Ginibre). 



feature. 



L 





Ck{X{g)) = Tk{^A9.9.---.9))- 
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Figure 3: A rotary flow: a cumulant term that does not vanish 

Looking to take advantage of the radial symmetry, we fix a simple test function C{z), a 
product of z and a radially symmetric weight function. Then taking the monomials fj to be 
powers of ( and its conjugate, 

m = (C(^))"^(cRF, 

we expand the corresponding integral to find that 

n „ m 

"^mifl, ■■■Jm)= / n Ci {Zj ) Ci {ZjY^ XjCj {Zj Z^+i) dfi{zi) ■ ■ ■ dfi{Zm). 

qi,...,qm=0 j = l 

Each term in the above sum can be thought of as a flow network; the vertices are 1 . . . m, 
and the directed edges are (j, j + 1), j = 1 . . .m (see Figure Ej). The exponent of Zj in the 
term is aj + qj, and the exponent of Zj is f3j + Qj-i, these have to agree or else the integral 
over Zj will vanish (granted by the rotation invariance of the reference measure). In flow 
network terminology, the node law has to hold (total inflow=total outflow), see figure El 
This implies that setting 

7i = /^j-ttj, j = l,...,m, 
Vj = 7i + --- + 7j, 

and also £ = g^, it is required that 

q. = i + rij, j = l,...,m, (9) 

ai + . . . + am = Pi + . . . + Pm, 

for the integral not to vanish. We refer to 7, /3, 77, q satisfying these conditions as a rotary 
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flow. Denoting 

M{q) = J^\z\''df^{z) = c^l^, 

Miq,K) = [ \zmz)rdfi{z), 
Jc 

we find the following master formula, 

^ — Vmin 3 1 

for the basic cumulant term. 

In the flow language, the conditions © say that for each edge (j, j + 1) in the cycle we 
get a term 

= XgM{2qjy^ = Xe+^M{2i + 2r],)-\ 

For the vertex j the exponent of the \zj\ factors coming from the kernels K{-, ■) equals the 
total flow within the cycle in or out of j: 

qj-i + qj = 2i + r]j_i + r]j = 2i + 2T]j - jj. 

The total exponent of |C(^j)| is the total external flow in or out of j, which equals aj + (3j. 

5 Combinatorial identities related to rotary flows 

The main goal of this section is to establish general conditions under which the combinatorial 
sums over rotary flows arising in connection with the limits of joint cumulants vanish. 

Let be a vector space; we will use = M unless otherwise specified. For a function 

a : [k] ^ [m], and a eV^ define aa G V™' by 

i:a{i)=j 

Fix a function '■ V"^ —>■ K for each 1 < m < k. Define 

(7:[m]—»[k] 

with the sum over all surjections, as well as 

Afc(v2, a) = ^k,m{<^m, a)- 

m 

m=l 
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The connection to the cumulant expressions above should be clear. By a shght abuse of 
notation, we will use the shorthand — for the collection of functions </?^ = af + ... + cc^. 
Also let 

Afe,m(l) = |{a : H ^ [A;]}| = Afc,„(x, (1,0,..., 0)), 
and introduce the shorthand Afc(l) = Afc(a;, (1, 0, ... , 0)). 
Lemma 5. For linear </? it holds that, 



Ak{x,a) 



"1, ifk = l, 
0, ifk> 2. 



Proof. The case A; = 1 is clear, so assume k > 2. Each (pk{oi) — cci + . . . + ctfc, which we 
denote by s. We expand each A^^^ according to whether a{k) has a unique pre-image. This 
produces 

— l,m— 1 

(1) + smAfe_i,^(l), 

and so the expression for Afc(x, a) is just a telescoping sum. □ 
Lemma 6. In the quadratic case we have 

aj, ifk = l, 
Kk{x^,a) = { 2aia2, if k = 2, 
0, ifk> 3. 

Proof. Again assume k > 2. Let a~ — {ai, . . . , CKk-i), and let s~ = cti + • • • + ak-i- We 
expand Afe,m in terms of a^, first considering the case where a{k) has a unique pre-image k. 
Summing over values of a{k) we find 

k {Ak-i,m-i{x^ : a") + alAk-i,m-i{^)) ■ 

Now take the case when a{k) has more than one pre-image. Expanding each quadratic term 
yields 

a:[k]-»[m] 

a (i) — (T(k) 

Performing the sum in the last term over the values of a{k) produces 2akS~ Ak-i^mi^) ■ It 
follows, 

Ak,m{x'^,oi) = mAk-i,m{x'^ , oi~) + mAk-i,m-i{x'^ , a') 
+ malAk-i,m{i) + malAk-i,m~i{i) 

+ 2s"Q!fcAfc_i,^(l), 
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and so, when we sum in m to compute A^, the first two rows telescope leaving, 

Afc(x^a) = 2s"afcAfc_i(l). 
This equals 2aia2 for k = 2, and for A; > 3 by Lemma El □ 



Lemma 7. Assume that each (fm, 1 < m < k is a (not necessarily symmetric) quadratic 
polynomial in the Ui. Let bm i = 0,1,2 be the sum of the coefficients of the a* terms in ipmi 
let bm^^ denote the sum of the coefficients of ajaj/ for j ^ j' . If 



6W, ^/A = (0), 

mb^'\ if\ = (1) or A = (2), 

m(m - if\ = {l,l). 



for some b^ and all m, then 



5(0) + + 5(2)^2^ z/A; = l, 

a) = { 2(6(2) _ 5{i,i))a,a2, if k = 2, 
0, tfk>3. 



Proof. Denote by 



- y 



nl 

a-.\n\^\n\ 



the symmetrized version of (fm- By definition A/^^rnifm, Ci) is symmetric in the a^. So, for 
a : [k] ^ [m] and y = aa, if we set s„ = yi + . . . + ym, then s^- = s = ai + . . . + ak does not 
depend on cr or m. Thus, ViVj = — YlT=i vh ^"^^ further 



m 

2 



i=l i=l 

From here we also have 

= ^b(o) ^ 5{i,i),2)^^^^(i) ^ 6(i)A,,„(x, a) + (6(2) - 6(i-i))A,,„(x2, a), 

and the claim follows from Lemmas El and El □ 

Lemma 8. Let now V = R^, and assume that each '■ R, is a (not necessarily 

symmetric) quadratic polynomial in {ai, Let 6^, 6~ denote the sum of the coefficients 
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of the aijSi and ai(3j, i j terms in ipm, respectively. Assume as well that all other coefficients 
vanish, and that there exists and h~ so that 

for all m. Then it holds 

h^ai(3i, if k = 1, 

Ak{v,(^)={ (6+-r)(ai/32 + a2/3i), ^f k = 2, 

if,k>3. 

Proof. By Lemma [7| we may assume that each is symmetric. Setting s^g = /3i + . . . + Pk, 
Sq, = «! + ... + ak, and for a : [m] -» [k] setting y = aa, w = aP, the arguments of Lemma 
[3 give 

m 

LPm{{yi,Wi), {ym.Wm)) = b'SaSfS + (6+ - 6") 

i=l 

We use polarization {^qjiWi = [yi + — {yi — WiY) for the second term and find that, 

The proof is finished by invoking Lemmas El and IHl □ 

Lemma 9 (Spohn, Soshnikov). Let (pmi.li, ■ ■ ■ , 7m) = max(r7i, . . . , rjm), where rjj = 71 + 
. . . + 7j . If ai + . . . + ak = 0, then 



Aki<^,a) 



— |ai| = — |a2|, if k = 2, 
0, ifk^2. 



A result of this type was first discussed by Spohn in |2I] (though see [201 as well) which 
focusses on the determinantal process on the line with the sine kernel SlsMizlll. it also 
appears as the 'Main Combinatorial Lemma' in the the paper of Soshnikov, [23j, where it is 
used to track fluctuations of the eigenvalues of the Unitary group U{n). 

In our result. Lemma IHl figures into the boundary component of the limiting noise, and 
so may be thought of as the root of the connection to U{n). 

We finish with the note that in the formulation of |2S1, the rjj ='-/[ + .. .+7j_|. That is, the 
sum only goes up to index j — 1). This is equivalent to the above upon setting 7^' = —•jk+i-i, 
so that r]'j = 7^ + . . . + 7^._^ = -(7^ + . . . + 'jk+i-ij^i)) = Vk+i-j and max(?7^, ...,r]'m) = 
max(?7i, ...,r]m). 
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6 Polynomial statistics for the Ginibre ensemble 



We now apply the results of the previous two sections to the Ginibre ensemble. That is, in 



we put L = n < oo, ^{z) = ^n{z) = -e "'^1 , Ck = jt, = 1, and then take n oo. 



Further, while the possibility of using weighted polynomials may be important in general, in 
the case of Ginibre it suffices to work with the naked monomials = z. 
With these specifications. 



and plainly 



M{2i + 2r]j + aj - (3j, aj + (3j) = M{2i + 2rij + 2aj), 
M{2q + 2K) n'i{q + K)\ 



M(2k) 

Substituting into (fTUj) . we have 



q+K 



n ''{q+l)---{q + k). 



n—l—r}max m 



, . . . , 



^ — Vfnin j — 1 



which for large n reads, 

n' 



1 " 

- E 

where again s = ^ = ^ Next, use the fact that 



+ 1 



+ o{e 



•s-2\ 



E 



£ = const. 

valid for any s > 1, to put the the asymptotics of (|TT|l into the form: 



$m(/l,...,/m) 



n 



s + 1 



- (1 + ^n 



1 1 



1 1 \ ^ 



+ 0{n-^). (12) 



Here fi = fe{z, z) is shorthand for z'^^z^'-. 

For fixed s, the 0(1) term in (jl2p is a function of the and which we denote ip„ 
Summing the above into the basic cumulant term, we find that 



T,($,(/i,...,/,)) 



n 



s + 1 



A,(l) + Afc(<^,(a,/3)) + 0(n-^) 



(13) 



and are in position to use the combinatorial tools developed in Section 

First of all, the coefficient of n in ()13|) vanishes for A; > 2 by Lemma El For the constant 
order term, examining iprn shows it to be a linear combination of rjmaxi a quadratic polynomial 
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in ttj and a polynomial which is itself a linear combination of the ai(3j. Further, the sums of 
the coefficients in the polynomial part are as follows: 



1 


ai 




aittj^i ^ j 




ail3j,i ^ j 


-1/2 


m/{2s) 


-m/(2s) 


—m{'m — l)/(2s) 


m/ s 


m{m — 1)/ (2s) 



Now combined, Lemmas [71 |H1 and 121 imply that Ak{ip, (a,/3)), and so the full cumulant, 
vanishes in the limit for k > 2. In summary we have: 

Theorem 10. For any real-valued polynomial P{z,z) = Cai^^p^z"''' z^'' , the Ginibre statistic 
Xn{P) — E[X„(P)] converges in distribution to a mean-zero Gaussian as n (yo. 

As for the limiting variance, choose a pair of monomials z°''^z^^ and z°'^z^'^ with ai + ^2 = 
Pi -\- P2, and note that our cumulant asymptotics enter the picture as in 



lim Cov(2"l^^^^^2^"2^ 



lim f$i(2°i+"2^^i+/32) _ ^^(z'^^z^K 

n— >oo \ 

max((J, pi — ai) H . 



(14) 



On the left, z°'^z^'^ and z^^z'^^ stand in for the full linear statistics, and 02, (32 switch rolls from 
left to right since covariance is conjugate linear in the second argument. The limit itself is 
read off from (|T2|l : why the corresponding covariance vanishes in the limit if ai + a2 7^ /5i + /?2 
should also be clear. 
Next observe that 



a2Pi a2Pi 



and 



- [ d{z''^z^^)d{zP^z^Ad'^z -- 



J2 k {z"'z'^')''{k) (z/52z"2)^(A;) = max(0, ai - (3i] 

k>0 



(15) 



(16) 



A little algebra will show that the sum of ()15|) and ()16p equals the final expression in p4|) . 
By linearity we may conclude: with any real- valued polynomials / and g in z and z, 



lim Cov{f,g) 



The general covariance/variance asymptotics occupy the next section 



7 Concentration 

This section is devoted to the following estimate. 
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Theorem 11. For linear statistics and Xn{g) in the Ginibre ensemble, 



lim Cov(X„(/),X„(^?)) = - / dfiz)dgiz)d'z + y2kfik)g{k). 



(17) 



A:>0 



as long as f and g possess continuous partial derivatives in a neighborhood of U and are 
otherwise bounded as in \f{z)\ V \g{z) \ < Ce^'^L 

Granting (fT7|) . the central limit theorem extends immediately from polynomials of the 
form 

afe+/3fe<A/ 

to general test functions f{z) satisfying the growth and regularity conditions assumed in the 
statement. 

First, by the Stone- Weirstrauss theorem, we can find a polynomial P'{z, z) such that 
WP'i z, z) — (i/||ioo('|2|<;^^5) < e for whatever e > 0. Certainly the L°°-norm (of the derivatives) 
on the larger disk controls both the H^(JJ) and H^^'^^dV) norms, and there exists a sequence 
of polynomials Pm{z, z) such that 



f-Pn 



0, 



(18) 



Hi{u)nHi/2(u) 

as the degree m — > oo. At each step, Pm is just the anti-derivative (in z) of the polynomial 
approximating df. 

Next, denote the centered statistic by 

n 



Z^{f) = X„(/) - - / f{z)d'z. 
With / as in the above result, a separate (and easy) calculation shows that 



n 



E[X„(/)] - - / f{z)d'z 



oil] 



and it follows that the family {Z„(/)} is tight as n ^ oo. Denote a given subsequential 
limit by Zoo{f)- For polynomial test functions, we have proved that Zn{Pm) converges 
in distribution to a mean zero complex Gaussian, Zoo{Pm)- Now pick a smooth bounded 
: C ^ [0, oo) with || V0||oo < 1, and note that. 



E0(Zoo(/))-E0(Zoo(P,n)) 



< lim sup E 



0(Z„(/)) - 0(Z„(P„)) 

< limsupE[|Z„(/-Pj|Y/', 



where things have been fixed so the right hand side tends to zero as m — > oo. This appraisal 
is independent of the special subsequence chosen in the definition of Zao{f)- Hence, by 
appealing once more to Theorem ^2 and ()18|) the large m limit of Z^{Pfn) identifies the 
limit distribution of Zn{f) unambiguously and completes the proof of Theorem [TJ 
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Verification of the limiting variance 

The opening move in the proof of Theorem ^2 is the covariance formula, 



Cov(X„(/),X„(^) = f{z)g{z)Kn{z,z)di2n{z) 



c 




f{z)g{w)\Kn{z,w)\ diJn{z)djj,n{w), 

vahd for any determinant al point process. Then, each appearance of the test functions / 
and g is expanded by way of the dbar representation. That is the fact that, for any / once 
continuously differentiable in a domain D C C and continuous up to the boundary 9D, it 
holds, 

;(C).-1/M4,.., 1 / /M,.. (19) 



see for example 

The formula ()19|) allows one to pass from the general variance asymptotics to one for 
the particular functional z [z — ■)~^. The decay properties of the Ginibre ensemble 
allow further simplifications. In particular, with f{z) once differentiable in say \z\ < 1 + e, 
decompose as in f{z) = f{z)ip{z) + f{z){l — ip{z)) for a smooth ip, equal to 1 on U and 
vanishing for 1^1 > 1 + e/2. Next, for any e > and c > 0, 



/ e''\'\'Kn{z,z)diiN{z) = / e^l^I'lzpV^I^I'd^z (20) 

J\z\>l+e J\z\>l+e 



n roo 

< nx / r^-^e-^"-'^)^ dr < C"r^e-""'/^ 

(n-1)! J(i+,)2 

by a simple Laplace estimate. From the growth assumption on / and noting \K{z,w)\'^ < 
K{z, z)K{w, w), it follows that the variance of Xn{f) agrees with that of Xnif'ip) up to o(l) 
errors. The same pertains to the covariance, and it suffices to take the test functions com- 
pactly supported from the start, ignoring the boundary integral in the dbar representation 
fll9|) . In short, the covariance may be read off from the n ^ oo behavior of 



7r2 



df{u)dg{r]) X (21) 

\u\<l+e J\'r}\<l+e 



— ^ Kn{z,z)dfi{z) - [ [ ^ ^ \Kn{z,w)\'^ dfin{z)dfin{w) \ d^ud^T], 

ICZ-U Z-T] JcJcZ-U W-T] J 

in which f{z) and g{z) and their first partial derivatives are assumed to vanish along \z\ = 
1 + e. 
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Defining, 



4>nii^,r]) = I — — Kn{z,z)d^n{z) (22) 

Z — U Z — 7] 




Kn{z,w)Kn{w,z)dfin{z)dfln{w) 



Z — V W — T] 

the proof of Theorem 1111 splits into two parts. For the interior (or H^) contribution we have: 
Lemma 12. For fixed v with < 1, 

in the sense of measures as n ^ oo. It follows that, 

lim \ [ [ df{u)dmM'^,n)d'ud'v=- [ \df{u)\'d\ 

for bounded and continuous df . 

For the boundary, (or H^f^) contribution, we prove separately: 

1 /2 

Lemma 13. With Hj^ (9U) corresponding to the usual sum extended only over positive 
indices, it holds that 

lim [ df{i.)W¥)M^,v)d'ud'v=\\f\\l.^.,,^y 

whenever f is continuously differentiahle and vanishes on\z\ = 1 + e. 

Note that we have set f = g in the above two statements. This has only been done for 
the sake of slightly more compact expressions. No generality is lost: covariances can always 
be recovered from variances. 

Proof of Lemma\T^ We begin by showing that, 

lim / 4>n{^,v) d^V = 71", (23) 
for < 1 and B any disk about the origin containing u. A residue calculation gives, 

d^rj J TTz for \z\ < b 
\<bZ-V~ \ ^ for \z\>b ' 
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and so, denoting the radius of B by 6 > \u\ = a, 

I (l)n{i^,f])(fr] = nj Kn{z,z)dfin{z) + 7ib'^ -^Kn{z,z) dfin{z) (24) 

-'a<|z|<b ~'|2|>fe 
n—2 „ „ 

» n J\z\>a J\w\<b 



n-2 



J|2|>a J\w\>b 

k 

after expanding (z — z/)~ in series where \z\ > |z/|. Recall that Ck = = Now recombine 
(j24p in the form 

[ M^,v)d^V = vrc^-i / \zf''-^Uiin{z)-^ f (l-^)diin{z) (25) 

JB ^|2|>a ^|2|>6^ 1^1 ^ 

-ttV'qq+i / / 1 - —f^ ) dyu„(z) (i/i„(tu). 

^ ii2|<ay|«,|>6 ^ kr^ 

A simple computation shows that the first two terms of ()25|) tend to vr and zero respectively 
as n — i> oo. The final two terms of ()25|) may in turn be bounded by a constant multiple of 



n—2 „ „ 71—2 

/ / |^PVP^+2ci/x„(z)d^„(«;) < J2p(s^ < a')p(s^^, > b') (26) 

£=0 J|2l<a J|i«|>fe £=0 



where denotes a sum of {£ + 1) independent mean 1/n exponential random variables. 
Standard large deviation estimates explain why ^(5"" < a^) is exponentially small in n if 
i > n{a? + 5), with an estimate of the same type holding for P{S^ > 6^) and i < n{b'^ — 5). 
Since there is a gap, i.e. a < b, each term in the above sum is of order e"'"'" for a 7 positive 
with (6 — a)^. Since there are only n terms, the verification of is complete. 

It remains to show that we have decay away from rj = u. Given 5 > 0, we claim that 



lim 

n— ►oo 



\\u—rj\>S 



d^r] = 0, 



with any A supported in |z| < 1 + £. 

First fix 77 as well as u, and take a = {ul < \r]\ = b without any loss. Two types of 
expressions emerge from performing the integrations over z and w in the definition of 0„(z/, fj), 
an inner term corresponding to \z\ < a and \w\ < b, and an outer term corresponding to 
\z\ > a and \w\ > b: 

0„(z/, r/) = (f)o{iy, fj) + (/)e(z/, r/). 
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The integration over the mixed regions < a, |w| > 6, for example) vanishes by orthogo- 
nahty. The inner term reads 



n— 1 ^ n—l „ 



|2(m+€) 



d\ln{z) 



(27) 



n-l 



£=n m=0 -^I^K" 



2(m+^) 



(i/i„(2;) 



n— 1 ^ n—l—\ 
^=0 ^ '-^ m=0 



The result of the outer integration has the same shape, up to the obvious inversions. For 
example, the analog of the first term on the right of (P7|) is 0^^'' = X]"=i YTmln-t i'^vY 
Cmi|2|>b \^\'^^'^~^~^^dfiniz). We sketch the proof of the decay of 0o. The estimates for 0e 
are much the same, and both cases are similar to considerations immediately above. 
Integrating by parts produces the bounds, 

1 



^ a^^+^ e~"°' < / Izl^^du (z)< 



^ 2k+2 ^-na^ 



(28) 



for k > n. Recall here that a < 1. If we further assume that b > a, there is a constant C 
such that 



^-1 „n-a-m) 
-na^ ^ 

{n-(£-m))\ 

£=0 m=0 ^ ^ 



2{n-(£-m)) 



a/6)' 



fc=0 



n—k 



k\ 



(after changing variables and the order of summation) and 

oo n—l 



/ ')\h 

na ) 



k=0 



k\ 



< c- 



[a/by 



1 - a/b)' 



The last two displays clearly tend to zero as ^ oo, providing the advertised L^-decay over 
sets supported away from = |?7|. To finish it is enough that both (j)'o\i',ff) and (j)i'\i',ff) 
remain bounded along that part of the circle \v\ = \ri\ where u/t] = e*^ and say \6\ > 6/2. 
This is a simple exercise taking advantage of the oscillations introduced in summing powers 
of e*^ 
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As for (0\v^r]), the estimate ()28p does not have the same effect as m + £ < n in each 
appearance of / \z\'^'^'^^^^ djjin{z) . Instead, one works along the hnes of pUj) . Again, first take 



a <h and write, 

n~l n-l-l 



=1 m=0 

.2 



Focussing on the sum restricted to £ + m < na , that object is bounded by 



no? na?—(. na^ na^—£ , ^ 

£=0 m=0 e=0 m=0 



r,2 



i=0 

having once more used the standard large deviations estimate for exponential random vari- 
ables. The sums over m + i > nb"^ and na? < m + i < nb"^ are handled by the same 
procedure. When a = b but u/rj is kept away from 1, the key observation are that 
^i^m+e ^ (^'^)PiSm+e ^ '^^) has a good decay away from \m + i — na^\ < C\fn. By 
using this in conjunction with the oscillations from [v/r])^^, the boundedness of (f)o^\i',fi) 
along {|z/| = |?7|} fl {{u — ri\ > 6} will follow. □ 

Proof of LemmaWA Start with the point-wise limit of 0„(z/, f/) for which we fix z/ 7^ in 
the annulus 1 < | ■ | < 1 + e. For the interior estimate, it was be convenient to carry 
the integration in z and w in the definition of 0„ over all of C, see fl22j) . However, by the 
estimate (j2Uj) . one may cut down the integral to | ■ | < 1 + 5 for any 5 > at the expense of 
an exponentially small error. We therefore consider the n 00 limit of 

4>n{v,f]) = I 1 Kn{z, Z)dfln{z) 

J\z\<l+5 {Z-J^){Z-V) 



|,|«,|<i+5 iz-iy){w-r]) 

in which 6 is chosen so that 1 + 5 is less than either or \ri\. 

Next, each appearance of (• — rj)^^ and (• — u)^^ in (pn is expanded in series, and the 
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integrals are performed term- wise to find that, 

c-n J\z\<l+5 K^V) 



n-l ^ ^ , 



/ , . \Kn{z,w)\'^d^n{z)dfin{w) 

^ J \z\<i+s J \w\<i+s \y^) 

(-_n ^'^Vj J\z\<l+5 



e=0 m=0 -'|2|<l+<5 

In the first equahty, only the diagonal terms survive in the expansion an account of orthog- 
onality. Now note that, for any a > 0, 

in which > for 5 > 0, and 



l-Ck f Izl^'di^niz) = 0(e-"(^-")') 

J\z\<l+5 



'\z\<l+5 

for any A; < (1 -|- a)n with a < 5. It follows that, 

\an\ n—1 



£=1 ^ m=n-£ 



with a obeying both the above constraints. Finally, since 

n-l e-1 



m=n-^ '^"*+^ m=0 



we have that 



lim (/>„(i/,77) = V 

r). — i-r\~\ • ■* 



The limit holds uniformly on{|i/| > l-|-(5u|?7| > 1-1-5} for any S > 0, and a dominated 
convergence argument produces 

lim 1 / df{u)dmM'^: V) d'ud'ri = E d - / ^ ci^r^ 



24 



It remains to realize that since d[f{r])ri ^{rj — uj)] = f] ^{f] — uj){df){ri) for r] away from zero 
and any cu, the dbar formula reads, 

= / ^TT (^V + TT -ITT ^V + TT -ITT 

^ Ji<\v\<i+e 27r r^^+i 27r J|^^|=i+e 

The third term vanishes by assumption, and the second term is precisely V— 1 times the 
£-the Fourier coefficient of the boundary data of /. □ 

8 Universality for analytic functionals 

For each positive integer n, pick a rotation invariant measure, dfj,n{z) = fin{\z\)(fz, and set 

n—l „ 

Kn{z,w) = ^Cn/{zidY, where c~^^ = M{n,2£) = / \z\'^'^dnn{z); 

n—n C 



the measure /i„ itself need not depend explicitly on n. Now impose the following moment 
condition: for all integer m there is a p > such that 

M(n, 2n + 2m) /M{n, 2n) p^"", (29) 

as n — > oo. Of course, by a scaling, it may be assumed that p = 1. 

For any determinantal point process {K, p) in C with rotation invariant p, it is the 
case that the moduli \z2\, ■ ■ ■ , \zn\} as a set have the same distribution as {Ri, . . . , Rk} 
where Ri are independent with P{Rk G dr) = 2TiCkr'^^^^ jl{r)dr for each k (see ^T] or 
Jl]). Therefore, ()29p has the interpretation that the stochastically largest random variable 
Rn satisfies EnR"^ —>■ p, for every m G Z. Equivalently, for every fixed £ > we have 
'EnRl-e P- 

Being the main example at hand, the reader will be happy to check that Ginibre satisfies 
everything asked for: the relevant computation is lim„^oo ^n"r(n) = 1- A second example of 
interest is the truncated Bergman ensemble. Here one begins with p^ = the uniform measure 
on U, producing the kernel 

n-l 

Kn{z,w) = ^(£+ 

e=o 

of orthonormal polynomials on the disk. This model arises naturally in the following way. 
Consider the random polynomial 2" + Y12=o ^kZ^ with independent coefficients drawn uni- 
formly from the large disk RU. If we condition the roots to lie entirely within the unit disk, 
the R ^ 00 limit of the resulting point process is the truncated Bergman ensemble. This 
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observation may be gleaned from Hammersely PU]; but see also ^7] for why the adjective 
"truncated" is used. 

Consider now a linear statistic in any such ensemble which is polynomial in z alone, or 
of the form YTi=\Vm{zi) = YlT=i^k(Yli=i4y We have the following CLT. 

Theorem 14. Let zi, . . . ,Zn be drawn from a determinantal process {Kn, fJ^n) (is above for 
which holds with p = 1. Denote pj{z^"') = zl + ■ ■ ■ + z^. Take a = (ai, . . . , a^) and 
b = {bi, . . . ,bk) with aj, bj G {0, 1, ... } and bounded independently of n. Then, 

k 



lim E 



n 



■ k 

l[pj{z^T'^;(^) 

.i=i 



E 



(30) 



for Zi, . . . , Zk independent standard complex Gaussian random variables. 



Having identified the limiting moments, this implies Theorem 3. On a case by case basis. 
Theorem IT^ and a variance estimate yields a more complete picture. 

Corollary 15. Let f{z) be analytic in a neighborhood of \z\ < 1 and the points Zi,. . . ,Zn 
be drawn either from the Ginibre or truncated Bergman ensemble. Then, as n ^ oo, 
S"=i /(-2^) ~ ^f{^) converges in distribution to a mean-zero complex Normal with variance 

Qualitatively, the Ginibre and truncated Bergman ensemble have marked differences. 
While the Ginibre eigenvalues fill the disk as n — oo, the truncated Bergman points con- 
centrate near \z\ = 1 {i.e., ^Kn{z,z) tends weakly to 6\z\=i). Once again, the analytic CLT 
only "sees the boundary" . 

The proof is largely inspired by the ideas of Diaconis- Evans ^ where the analogous 
moment formula is established for the eigenvalues of the Haar distributed unitary group. 
Those points yield yet another example in the above class: there dfin{z) = 6\z\=i and p = 1. 
Also in that case, jB] shows the equality (jHUjl to hold for finite n as soon as n > (X]^=ii%) ^ 
iYl^=i j^j)- The basic observation is that the integrand on the left hand side of (plj) is 
comprised of symmetric polynomials in the points zi,. . . ,Zn', one would like to to expand 
this object in a convenient basis. 

Let An be the vector space of symmetric polynomials in the variables z®" = zi, ...,Zk of 
degree at most n. Let A„ denote the set of partitions of integers at most n. Given a partition 
A = (Ai > A2 > • ■ ■ > ...) G A„, bring in the corresponding Schur function, 

det[zt^^-- 



Sxiz^ = \ .n''=\ (31) 

kl=l 



det izl 
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and recall the following well known facts. 



Theorem 16. ('WJ, Chapter 1) The Schur functions sx{z®'^)^ \ G A„ form a basis for A^- 
Consider the inner product (-,■) that makes them an orthonormal basis. For fn{z®"') = 
n^=iPfc'°(-2®") (ind (7„(-2®") = Y['j=iP''ki^^^) ^ ^n, which are products of simple power sum 
functions, it holds that 

n 

{fn,9n)=Sab{[[pa,\). (32) 

The inner products are compatible as n varies in the sense that if m > n, and X E An we 
have {fn,sx{z'^n) = {fm,sx{z'^n)- 

Lemma 17. For Schur functions in the points of {K, fin) o-s above, 



n-l 



£=0 



M(n,2(n + A£-^)) 
M(n,2(n-£)) ' 



(33) 



Proof. The denominator in (jHTj) is a Vandermonde determinant, and thus the interaction 
term (see (^) in the integral of question is cancelled: 



sa(z®")s^(z®-) 



det(z^^"^" ^) det(2:^''^" ^^dfin{zi) . . . dfin{z„ 



The normalizer here is Zn = n\YYe^iM{n,2i). Now, expanding each determinant on the 
right hand side, the generic term we get is a constant multiple of 



n „ 

1,-1 Jc 



k=l 

with a permutation a. Since A and vr are monotone, this will vanish by orthogonality of z 
and z when A 7^ vr. When A = vr, there is a non-zero contribution from the diagonal, a = id. 
To conclude, we compute 



E^ 



Sxiz^n Sxiz^"^)] =^ll ljz\'^''^--'Ufiniz) 



k=0 



but this is exactly (j33p . 



□ 



Proof of Theorem\T^ Let fn = 11^=1^^(2®")"^ and ^„ = 11^=1 Pi(-^®")^^ and let no be at 
least of the degrees of / and g (this is independent of n). Now expand with respect 

to the basis given by the Schur functions, and compute the expectation term by term. This 
expansion is finite, and the coefficients (/„, sa(-2®")), sa(-2®")) do not depend on n for 
n > no by Theorem ITBl 



27 



Since ()33|) has a bounded number of factors that are not 1, our assumption and Lemma 
IT7I imphes that for A, fixed we have E„[saS;/] — > Thus we have E„[/„^„] {fnoiQuo)- 
The claim now follows by ()32|) and the moment formula for complex Gaussians. □ 



Proof of Corollary{T^ For analytic f{z), /(O) = ^ j^f{z)(Pz = ^ Jg^f{z)dz, which ex- 
plains the centralizer. While for the Ginibre ensemble, we could simply quote the result of 
Section [7| the analyticity allows for a simpler approach amenable to more general ensembles 
subject perhaps to additional conditions in the spirit of (|^. Preferring to be concrete (and 
brief), we restrict ourselves to the truncated Bergman case. 

Rather than using the dbar representation, the variance in the analytic case may be 
computed by the more familiar Cauchy integral formula: with the reference measure now 
uniform on the unit disk, 
1 



47r2 



f{u)f{7]) I ^:,JL ^, dh- [ [ J^^^i^lp^d^zd'w dvdn. 



V 



{z - u){z -f]y ~ j^j^{z-v){z-fi)' 

Here is a circle of radius 1 + 5 > 1 (about the origin) within the region of analyticity of 
/. What to do next is plain: 

K„{z,z) 2 1 ^ m + 1 



-d'z 



oo ^ n~l 



and 




n~l , n-£-l 



-d zd w 



after expanding both the kernel and the functions {z — ■)^^, {w — ■)~^ and integrating term- 
wise. With \i>T]\ > 1 there is enough control to pass the limit inside the summations and 
conclude that the variance tends to 

as n — s> oo. Here we have used the analyticity to pull the integral back to dU after the 
limit was performed. This last expression, and much of the method getting there, is now 
recognized from the Ginibre boundary case fLemma I13|). In particular, that it equals the 
advertised ^ \f'{z)\'^d'^z has already been explained. □ 
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